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Abstract 

Mixed compression inlets are common among supersonic propulsion systems. However they are 
susceptible to total pressure losses due to shock/boundary-layer interactions (SBLFs). Because of their 
importance, a workshop was held at the 48th American Institute of Aeronautics and Astronautics (AIAA) 
Aerospace Sciences Meeting in 2010 to gauge current computational fluid dynamics (CFD) tools abilities 
to predict SBLI’s. One conclusion from the workshop was that the CFD consistently failed to agree with 
the experimental data. This thesis presents additional CFD and numerical analyses that were performed 
on one of the configurations presented at the workshop. 

The additional analyses focused on the University of Michigan’s Mach 2.75 Glass Tunnel with a semi- 
spanning 7.75 degree wedge while exploring key physics pertinent to modeling SBLI’s. These include 
thermodynamic and viscous boundary conditions as well as turbulence modeling. Most of the analyses 
were 3D CFD simulations using the OVERFLOW flow solver. However, a quasi-ID MATLAB code 
was developed to interface with the National Institute of Standards and Technology (NIST) Reference 
Fluid Thermodynamic and Transport Properties Database (REFPROP) code to explore perfect verses 
non-ideal air as this feature is not supported within OVERFLOW. Further, a grid resolution study was 
performed on the 3D 56 million grid point grid which was shown to be nearly grid independent. Because 
the experimental data was obtained via particle image velocimetry (PIV), a fundamental study pertaining 
to the effects of PIV on post-processing data was also explored. 

Results from the CFD simulations showed an improvement in agreement with experimental data with 
certain settings. This is especially true of the v velocity field within the streamwise data plane. Key 
contributions to the improvement include utilizing a laminar zone upstream of the wedge (the boundary- 
layer was considered transitional downstream of the nozzle throat) and the necessity of mimicking PIV 
particle lag for comparisons. It was also shown that the corner flow separations are highly sensitive to 
the turbulence model. However, the center flow region, where the experimental data was taken, was not 
as sensitive to the turbulence model. Results from the quasi-ID simulation showed that there was little 
difference between perfect and non-ideal air for the configuration presented. 
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1 Background 


1.1 Introduction 


Mixed compression inlets are often found in supersonic aircraft propulsion systems with a typical mixed 
compression inlet shown in Fig. 1. These inlets are susceptible to shock boundary-layer interactions (SBLI) 
with a typical SBLI shown in Fig. 2. As one can see, SBLI’s are not trivial in nature and beyond the 
situation shown in Fig. 2, are very three dimensional flows, shown in Fig. 3. SBLI’s are a concern for 
supersonic aircraft propulsion system designs as the separation induced by the shock interacting with the 
turbulent boundary-layer results in a total pressure loss that is not accounted for in inviscid theory. As mixed 
compression inlets become more at the forefront of aviation technology, it becomes crucial to understand 
SBLI’s and how they affect propulsion system performance. This can be achieved in several ways, including 
fundamental and flow control experiments as well as computational work, all of which will be discussed 
further in the next sections. 



Figure 1: A supersonic mixed compression inlet (used with permission from Pitt Ford and Babinsky [1]). 



Figure 2: Sketch of the oblique shock / boundary-layer interaction (used with permission from Touber and 
Sandham [2]). 
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Figure 3: Iso-surfaces of density gradients in the downstream direction around the wedge, colored by the 
derivative of the density gradients in the downstream direction (from Galbraith [3]). 

In order to determine how effectively computational fluid dynamics (CFD) tools can currently pre- 
dict SBLI’s, a workshop was held at the 48th American Institute of Aeronautics and Astronautics (AIAA) 
Aerospace Sciences Meeting featuring an array of CFD analyzes [14, 15, 16] with experimental data sets 
provided by the Institut Universitaire des Systemes Thermiques Industriels [17] and the University of Michi- 
gan (UM) [18]. The CFD results from the workshop are summarized in DeBonis et al. [19] with workshop 
conclusive remarks documented in works by Benek et al. [20, 21] as well as Hirsch [22]. One key conclusion 
of the workshop was that the CFD analyses failed to match experimental data. Because of the complex 
nature of SBLI’s, it is of great interest within the aero-propulsion and CFD communities to see if there are 
ways to improve CFD methods or more effectively use existing methods. This will enable better prediction 
of SBLFs as well as aid in future inlet design. 

Thus further CFD analyses were performed at the University of Cincinnati’s Gas Turbine Simulation 
Laboratory as a compliment to work done by Galbraith [3, 23] in hopes to better understand the compu- 
tational factors involved in calculating SBLFs as well as to explore alternatives to the error metric used 
in the workshop. One of the data sets provided to the workshop by UM was their Glass Tunnel with a 
semi-spanning 7.75 degree wedge at a freestream Mach number of 2.75. The UM experimental data was 
obtained using stereo particle image velocimetry (PIV) techniques [18]. The CFD analyses presented here 
focus on that particular UM case while exploring various key physics associated with SBLFs, including, 
but not limited to, heat transfer boundary conditions, geometry sensitivities, laminar verses turbulent flow 
assumptions, and turbulence modeling. Special attention was paid to the u and v velocity components, 
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particularly around the oblique shock off the leading edge of the wedge. This is because it was felt that the 
CFD and post-processing calculations from the workshop missed the peak u velocity as well as the location 
of the shock as defined by the v velocity profile. If the upstream profiles are not in agreement, it is fair to 
say that the solutions downstream should not agree. Aside from the CFD cases, a preliminary exploration 
into the effects on PIV obtained data is also presented. 

1.2 Experimental Research Literature Survey 

The workshop was not the first and certainly not the last time SBLI’s have been examined both experimen- 
tally and numerically. Experiments by Holden and Babinsky [4] attempted SBLI flow control via a series of 
streamwise grooves and bumps, shown in Figs. 4 and 5, with hopes of smearing the shock footprint. They 
showed that the flow fields inherent with SBLI’s are highly sensitive to geometry deviations, such as the slot 
and bump geometries, or for that matter any foreign debris or geometry imperfections. Another experiment 
by Holden and Babinsky [24] explored the use of vortex generators (VG’s) for SBLI flow control. Two types 
were tested: wedge-shaped, more commonly known as micro-ramps, and counter rotating vanes. It was 
found that the use of the VG’s greatly reduced the separation region of the SBLI interaction region, with 
the vane type VG’s going as far as to eliminate the separation region completely. Experiments by Pitt Ford 
and Babinsky [1] as well as by Lapsa [12] also explored the effects of micro-ramps on SBLI flow fields. Pitt 
Ford and Babinsky showed that micro-ramps located upstream of the interaction region can break up the 
separation bubble but not completely eliminate it while increasing downstream velocities, shown in Fig. 6. 
Lapsa showed that the use of inverse micro-ramps can decrease the displacement thickness and thus allow 
for less separation around the interaction region. 


a) Datum groove - 
4mm wide, 20mm deep 

b) 8mm wide groove - 
20mm deep 

c) Tapered nose 

d) Tapered tail 

e) Filledl 

f) Filled2 

g) Single baffle 

h) Two baffles 


Figure 4: Datum groove, 8mm wide groove and variations (from Holden and Babinsky [4]). 
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a) Wedge 


b) Ramp 


c) Bumpl 





Figure 5: Three-dimensional bumps (from Holden and Babinsky [4]). 



Figure 6: Velocity profiles 94 mm downstream of micro-ramp controlled SBLI, h = 3 mm (used with per- 
mission from Pitt Ford and Babinsky [1]). 

While the above mentioned experiments focused mostly on flow control, more fundamental experiments 
have been performed in order to understand the nature of SBLFs. An experiment performed by Babinsky 
et al. [25] explored the importance of corner flows in relation to SBLFs. They showed that the corner flow 
separations have a coupled effect with the separation bubble of the interaction region of SBLFs. Thus by 
decreasing the corner flow separation it was shown that the interaction region separation would be reduced 
and approach a 2D nature. This coupling nature was also shown in a prior experiment by Titchener et al. 
[26]. 

1.3 Computational Research Literature Search 

Over the years there have been plenty of numerical simulations exploring all facets of SBLFs: from flow 
control to numerical modeling efforts. For flow control, Ghosh et al. [5] investigated the use of aeroelastic 
mesoflaps via two and three dimensional simulations, shown in Fig. 7. The simulations utilized both 
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Reynolds- Averaged Navier-Stokes (RANS) and Large-Eddy /Reynolds- Averaged Navier-Stokes (LES/RANS) 
means for solving the flow equations. The simulations showed that the use of mesoflaps for flow control 
resulted in a slightly larger interaction region compared to the non-controlled case. 



x (m) 



Figure 7: Mach number contours of the converged mesoflap array simulations, (used with permission from 
Ghosh et al. [5]). 


Morgan et al. [6] performed LES simulations with a flat plate geometry to explore SBLEs. An instanta- 
neous as well as time-averaged snapshot of the skin friction coefficient within the interaction region is shown 
in Fig. 8. It can be seen that the time-averaged snapshot greatly smooths out the chaotic nature of the 
SBLI, including the separation region. In fact, the instantaneous snapshot shows that the separation region 
consists of many separation bubbles and not just a single separation zone as shown in the time- averaged snap 
shot. This is important as steady RANS simulations only have the capability to reproduce time- averaged 
behavior unlike LES instantaneous solutions. The comparison between LES and RANS solutions were fur- 
ther explored by Medic et al. [7] using a simple square duct with results shown in Fig. 9. It can be seen that 
the RANS solution over predicts the boundary layer growth at the duct corners relative to the LES solution. 
This over prediction can have a tremendous effect in SBLI simulations, especially when using a tunnel-like 
geometry, due to the coupling effects between the corner flows and interaction region separation discussed 
earlier. Other LES solution efforts have been conducted by Hunt and Nixon [27], based on experimental 
data obtained by Dolling and Murphy [28], as well as by Jamalamadaka et al. [29]. 
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Figure 8: Contours of the skin friction coefficient of a flat plate SBLI with time averaged (top) and instan- 
taneous (bottom) snapshots (from Morgan et al. [6]). 


Streamwise mean velocity 

RANS SST LES WALE 

! ! 



yD yD 

Figure 9: Comparison of RANS and LES results on a square duct (from Medic et al. [7]). 

Further, a series of simulations were performed by Knight et al. [8] exploring SBLFs. Five configurations 
were run, shown in Fig. 10, and utilized Direct Numerical Simulations (DNS), LES, and RANS. There were 
a multitude of lessons learned from these simulations, especially regarding the RANS simulations. With 
RANS, the case of the 3D single fin correctly predicted the secondary separation region with use of the 
Wilcox-Durbin model. However, the RANS simulation of the 3D double fin failed to accurately predict 
the surface heat transfer using the linear and weakly non-linear Wilcox-based models. This is common for 
RANS simulations of SBLI flows, as also shown by simulations conducted by Knight and Degrez [30]. Other 
overviews of SBLI simulations involving DNS, LES, and LES/RANS methods can be found in works by 
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Zheltovodov [31] and Edwards [32], the latter of which touches briefly on experiments exploring the effects 
of heat transfer on SBLI’s by Zheltovodov et al. [33, 34]. 



2-D Expansion compression corner. 





Figure 10: Configurations run with DNS, LES, and RANS (used with permission from Knight et al. [8]). 


Wu et al. [9, 35, 36, 37, 10, 38, 39] have performed a series of DNS simulations on a 24 degree compression 
ramp with experimental data provided by Bookey et al. [40]. Originally they showed that the DNS simu- 
lations under predicted the separation region within the SBLI interaction, as shown in Fig. 11. Numerical 
code bugs as well as the incoming conditions to the corner were eliminated as factors of the discrepancy 
due to the agreement between the DNS simulation and the experimental data upstream of the compression 
corner. In turn, they determined that numerical modeling of SBLEs is sensitive to the numerical dissipation 
and suggested the use of a limiter on the dissipation schemes. Rerunning the DNS simulations with the 
limiter showed great improvements in predicting the size of the separation region, shown in Fig. 12. 

Continuing with the compression corner theme, Edwards et al. [11] used LES/RANS to explore a 28 
degree compression corner flow interaction. Menter’s Shear Stress Transport (SST) [41] turbulence model 
was used for the RANS portion. Although the simulations captured the shock behavior well, it was noted 
that there was a great dependence on the shear stress transport limiter within the SST turbulence model, 
shown in Fig. 13. Further analyses of the SST and other k-omega based turbulence models as related to 
SBLTs have been performed by Georgiadis and Yoder [42] as well as Tan and Jin [43]. One draw back 
to LES/RANS method used was the need to calibrate the blending function constant per case. This was 
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eliminated by Giesking et al. [44, 45] by using an estimate of the outer-layer length scale based on the resolved 
turbulent kinetic energy, ensemble-averaged modeled turbulence kinetic energy, and ensemble-averaged and 
time-resolved turbulence frequencies. This estimated length scale, in conjunction with the inner-layer length 
scale, was then used to determine to blending function model constant. 





Figure 11: Flow pattern sketches for the experiment (left) and DNS (right) for a 24 degree compression 
corner (from Wu et al. [9]). 



Figure 12: Comparison of the size of the separation bubble with (new) and without (old) a dissipation limiter 
(from Wu and Martin [10]). 
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Figure 13: Skin friction distribution (RANS models) for a 28 degree compression corner at Mach 5 (used 
with permission from Edwards et al. [11]). 
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2 Geometry and Numerical Modeling 


2.1 Geometry and Mesh 

The UM Glass Tunnel [12] is a suck down tunnel shown in Figs. 14 and 15 with the oblique shock generating 
wedge shown in Fig. 16. Although the upstream converging-diverging nozzle is interchangeable, only runs 
with the Mach 2.75 nozzle are explored in this thesis. The wedge is centered about the center-span of the 
tunnel and the tunnel sits in a room controlled to a temperature of 295.7 K ±1 K [46]. At the freestream 
velocity of the tunnel, the static temperature is about 118 K. The top and bottom walls are made of aluminum 
while the side walls and bottom window are composed of glass. The test section was designed to be 2.25” x 
2.75” with a throat cross-section of 2.25” x 0.742”. However, the current “as installed” dimensions deviated 
from this with a test section and throat cross-section of 2.25” x 2.72” and 2.25” x 0.725”, respectively [46]. 
These measurements were taken well after the experimental data had been collected for the workshop and 
the tunnel has been taken apart and reassembled since then. As such, it only offers an approximation of 
what the dimensions might have been for those runs. 



Figure 14: Cut-away view of the University of Michigan Glass Tunnel (from Lapsa [12]). 
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SEEDER 


FLOW STRAIGHTENING 


NOZZLE & TEST SECTION 



Figure 15: Schematic of the University of Michigan Glass Tunnel with upstream flow straightener and seeder 
(from Lapsa [12]). Note units are in inches. 



Figure 16: 7.75 degree wedge geometry. Dimensions courtesy of Lapsa [12]. 

The tunnel, along with the wedge, was modeled by a modified version of a 3D over-set grid by Marshall 
Galbraith. The original grid, containing 53 million grid points divided into 15 zones, paid particular attention 
to the packing along the walls and around the oblique shock location. However, it was felt that the throat 
region could benefit from a more dense axial clustering. Thus an additional 50 axial points were inserted to 
define the throat contour, shown in Fig. 17. The final grid of 56 million grid points is shown in Fig. 18. Like 
the original grid, the final grid was packed very tightly to the wall such that = 0.25 at the first point off 
the wall, based on fully expanded tunnel conditions at Mach 2.75. 

The grid coordinate system uses a left-handed coordinate system non-dimensionalized by the tunnel 
height of 2.75 in. However, the coordinate system used for data comparison is consistent with the one used 
in the workshop, which was dimensionalized in mm with the origin at the strut leading edge, bottom wall, 
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and center-span. To convert from the grid coordinate system to the data comparison coordinate system: 


x = 25.4 (2.7bx gr id - 38.265) (1) 

y = (25.4 x 2.75) z gr id (2) 

* = 25.4 {2.7by grid - 1.125) (3) 

Note that the grid coordinate system has the z-coordinate and y-coordinate flipped relative to the data 
coordinate system. It should also be noted that the data coordinate system is slightly different than the 
one used by Lapsa [12], in which the axial origin was at the inviscid shock impingement location. Both are 
different than the coordinate system used by Eagle et al. [46] for the planned second SBLI Workshop, which 
is a left-handed coordinate system with the origin located at the wedge leading edge, bottom wall, right wall 
(looking downstream). 



Figure 17: Grid throat region packing with original (left) and final (right). 



Figure 18: Side view of the grid representing the tunnel and 7.75 degree wedge with center-span wedge and 
throat inserts shown. 
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2.2 Solver 


OVERFLOW [47, 48], a RANS flow solver for structured over-set grids, was chosen as the CFD solver for 
the SBLI analyses. Version 2.2e was used on a cluster of 20 Quad-Core Xeon X5570 processors (NASA 
Pleiades-Nehalem). Runtime for each analysis took approximately three days on the 80 processors. For 
solving the Navier-Stokes equations, spatial integration used the HLLC scheme [49] with the Koren limiter 
[50] to third-order accuracy while temporal integration used SSOR [51] to first-order accuracy. All analyses 
utilized the SST turbulence model with all zones considered turbulent, unless otherwise specified. 

OVERFLOW requires a set of reference conditions to non-dimensionalize the input parameters and set 
the freestream conditions. Because OVERFLOW uses the freestream conditions to initialize the flow field 
when not given an initial solution, Mach 2.75 conditions were chosen for the reference state with total 
quantities equaling the room statics due to the vacuum driven nature of the tunnel. Using ID perfect gas 
equations yields the following reference values: 


T t = 295.7 K (4) 

p t = 98kPa (5) 

M re f = 2.75 (6) 

L ref = 0.06985 m (2.75 in) (7) 

7 = 1.4 (8) 

R = 287 J/ (kg . K) (9) 

Too = T, (l + 1 = 117.69 K (10) 

Poo = Pt (l + - 2 1 M^ef ) = 3.898 kPa (11) 
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y/'yRToo = 217.46 m/s 


( 12 ) 


^oo — 


Voq = M re fd oo = 598.01 m/s 


(13) 


Moo 


1.716xl0 -5 


Poo = = 0.1154kg/ 

oo 

Too V 5 / 491.6 + 198. 6\ 
491.6 ) V ^ + 198.6 ) 


m 


3 


= 8.163 x 10 6 kg/(m.s) 


(14) 

(15) 


Re L = P°° V °° L ~f = 592, 876 (16) 

Moo 

Note the test section height was chosen as the reference length due to the grid being non-dimensionalized 
by this parameter and was calculated using Sutherland’s Law. Specifics on the direct OVERFLOW 
input variables and redimensionlizing scheme are outlined in reference [3]. 
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3 Case Overviews 


3.1 CFD Cases with OVERFLOW 

To explore the key physics mentioned earlier, seven cases were run as outlined in Table 1. A baseline case 
(denoted as Standard) was first run with various parameters set to reflect prior CFD analyzes. To explore 
heat transfer boundary conditions, a case was run such that the top and bottom walls (along with the wedge) 
were isothermal while all other surfaces (including the bottom window) remained adiabatic. The window was 
approximated from 45 mm to 140 mm axially and -16.575 mm to 14.425 mm spanwise [46]. The isothermal 
surfaces were set to the room temperature of 295.7 K. Although the actual temperature of the aluminum 
varies due to heat transfer, assuming it is isothermal at the room temperature is likely to be better than 
considering it to be adiabatic. 

To explore geometry sensitivities, a case was run with the grid modified to reflect the as-installed geometry 
with greatest measured tolerances accounted for. The grid was modified by raising the bottom wall to achieve 
the desired test-section height and redefining the nozzle curve to obtain the desired throat height. The 
nozzle curve modification was achieved by dividing the contour into two sections and scaling accordingly. 
The upstream section was defined from the trailing edge of the inlet straight section to the geometric throat 
while the downstream section was defined from the geometric throat to the leading edge of the tunnel straight 
section. The new upstream contour section was then defined as: 


Vnew ( ythroat,new Vthroat^old) 


X\ — X 
*^1 % throat 


' Hold, 


(17) 


where 


m = 0.5 + 1.5 ( Xtfcroa * * ) (18) 

\ % throat %1 ) 

Note the subscript 1 denotes the inlet trailing edge match point. Likewise, the downstream contour section 
was defined as: 


Vnew — ( ythroat,new Uthroat,old) 


x 2 - X 

% throat 


T Hold, 


(19) 


where 


m = 1 + 2 


( % throat % 


( 20 ) 


\ % throat %2 , 

Note the subscript 2 denotes the tunnel straight section leading edge match point. A sample modified nozzle 
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contour with an exaggerated throat area based on this scheme is shown in Fig. 19. 



Figure 19: Sample modified nozzle contour. Note the difference between the old and new throat heights is 
exaggerated compared to what was actually used. 

To explore laminar verses turbulent flow, two cases were run. A base case was run with laminar flow from 
the nozzle plenum inlet to the leading edge of the wedge in order to establish a “trip” location for the second 
case. In OVERFLOW, the effect of a boundary-layer trip was simulated by forcing upstream zones to be 
laminar. OVERFLOW handles laminar zones by zeroing out the production terms of the turbulence model in 
use [47]. At the desired trip location, or transition point, the production terms within the turbulence model 
are activated for all zones downstream of this point. As a result, eddy viscosity is calculated, and added to 
the laminar viscosity in the transport equations. Figure 20 shows the momentum thickness Reynolds number 
contour for the bottom wall (half-span) while Fig. 21 shows the corresponding values at the bottom wall 
center-span. The trip location was defined where the momentum thickness Reynolds number approximately 
equaled 400. This value is suggested by Abu-Ghannam and Shaw (AGS) [52] as the start of transition for 
turbulence intensities of 1.5%; however several values of momentum thickness Reynolds number could have 
been chosen to define the trip location. Using a momentum thickness Reynolds number of 400 yields a trip 
location at an axial location of -470 mm. Thus the second case (denoted as Trip) had laminar flow defined 
from the nozzle inlet up to that trip point. OVERFLOW uses a grid line to set transition, so the gridline 
located at the trip point on the bottom wall was used as the trip location for all four walls. From Fig. 21, it 
can be seen that the momentum thickness Reynolds number is approximately 260 at the nozzle throat. The 
turbulence intensity must be greater than 2.3 for this flow transition to occur upstream of the throat based 
on the flat zero pressure gradient AGS [52]. Considering this and the favorable pressure gradient situation 
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of the accelerating nozzle flow, which in turn tends to delay transition, it is likely that the flow transitions 
downstream of the throat. 



Chosen Trip Location 

Figure 20: Momentum thickness Reynolds number at the bottom wall (half-span) for the laminar case. 


Momentum Thickness Reynolds Number for the Bottom Wall, Center-Span 



X (mm) 


Figure 21: Momentum thickness Reynolds number at the bottom wall, center-span. 

A case was also run that combined the attributes of the isothermal, as-installed geometry, and trip cases 
to demonstrate the combined effects each parameter has on the flow field. In addition, a sweep of cases was 
run to explore sensitivities to turbulence quantities and turbulence modeling. Two of these cases explored 
sensitivities to the freestream turbulent kinetic energy and freestream turbulent viscosity by running the 
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Standard case with higher values of each. The Standard and Combined cases were also run with an array of 
varying turbulence models, of which a turbulence model case study was performed prior to running the SBLI 
configuration and discussed later in this thesis. It should be noted that one sensitivity already accounted for 
in the above cases is the sensitivity to total temperature. The total temperature used for the workshop CFD 
cases was 293 K, a value that was not experimentally measured at the time. Using ID perfect gas equations 
yields a velocity increase of 2.8 m/s (or 0.47% of a 600 m/s freestream velocity) for the 2.7 K increase in total 
temperature. Further, the ±1 K in the total temperature measurement would yield an additional fluctuation 
of ±1 m/s (or 0.17% of a 600 m/s freestream velocity). These are all relatively small numbers, however, 
additively they may have an important combined effect. 

Table 1: Summary of CFD Cases. 


Case 

A/A* 

Turbulence Model 

Heat Transfer Boundary 
Conditions 

Freestream 
Turbulent 
Kinetic Energy 
(m 2 /s 2 ) 

Normalized 

Freestream 

Turbulent 

Viscosity 

Standard 

3.7062 

SST 

All surfaces adiabatic 

3.576xl0 -1 

0.3 

Isothermal 

3.7062 

SST 

Top/bottom walls and 
wedge isothermal at 
T=295.7 K. All other 
surfaces (including 
bottom window) 
adiabatic. 

3.576xl0 -i 

0.3 

Modified 

Geometry 

3.7847 

SST 

All surfaces adiabatic 

3.576xl0 -1 

0.3 

Trip 

3.7062 

Laminar from nozzle 
inlet up till trip, 
SST downstream of 
trip. Trip location 
at x=-470 mm 

All surfaces adiabatic 

3.576xl0 -i 

0.3 

Combined 

3.7847 

Laminar from nozzle 
inlet up till trip, 
SST downstream of 
trip. Trip location 
at x=-470 mm 

Top/bottom walls and 
wedge isothermal at 
T=295.7 K. All other 
surfaces (including 
bottom window) 
adiabatic. 

3.576xl0 -1 

0.3 

TKE 

3.7062 

SST 

All surfaces adiabatic 

3. 576x1 O'-' 

0.3 

MUT 

3.7062 

SST 

All surfaces adiabatic 

3.576xl0 :i 

3.0 


3.2 PIV Exploration 

To better compare the CFD with the experimental data, the CFD solutions from the Standard and Combined 
cases were re-post-processed to explore particle lag that is associated with the PIV techniques used to acquire 
the experimental data. To account for the particle lag, a crude model adjusted the coordinates of each CFD 
data point using the following equations: 
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X = X ± l^T 


(21) 


y = y + v v t (22) 

Three time constants, r, were chosen to represent a 50%, 75%, and 100% total velocity reduction shown in 
Fig. 22. This yields time constants of 1.8 (short), 3.7 (medium), and 5.5 (long) ps, respectively. Because 
the model was only applied to the center-span and the spanwise velocity is nearly zero, only the axial and 
transverse directions were accounted for in the particle lag calculation. Note that this model only attempts 
to mimic the effect of particle lag and that it is not based on the actual physics under consideration. To 
complement the particle lag model, a window averaging scheme was initially explored using a similar method 
used by Garman, Visbal, and Orkwis [53, 54] to average out the values within the center spanwise plane. 
However, the PIV grids proved to be nearly as fine as the CFD grids for the method to be effective with 
the prescribed window of 0.24 x 0.24 mm [46]. Window averaging was also attempted by averaging three 
spanwise planes, which included the centerspan and ±0.75 mm to cover the 1.5 mm spanwise differance [12]. 
This was shown to have little effect; the results of which are not shown in this thesis. 



n (mm) 

Figure 22: Measured particle response through an oblique shock (from Lapsa [12]). The velocity component 
normal to the shock, u n , is normalized by the pre-shock (u n i) and post-shock {u n 2 ) velocities and shown as 
a function of the shock-normal direction, n. An exponential fit to the data reveals the particle relaxation 
time, rp = 5.5 ps. 
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4 Flat Plate 


4.1 Overview 

Before exploring turbulence model sensitivities with the complex SBLI cases, a fundamental study was 
performed to gauge the effectiveness of various turbulence models available. The 2D Zero Pressure Gradient 
Flat Plate Verification Case provided by the Turbulence Model Benchmarking Working Group [13] was 
chosen for this study. Boundary conditions along with the finest grid are shown in Fig. 23. The flow across 
the plate is at Mach 0.2 with a Reynolds number of 5 million per foot and a reference temperature of 540 
°R. The plate itself is modeled as an adiabatic solid wall. Because OVERFLOW can only handle 3D grids, 
the 210,000 point 2D grid was converted into a 3D grid by duplicating the grid in three depthwise planes, 
for a total of 630,000 grid points. The 2D boundary condition was then applied within OVERFLOW to 
handle the extra dimension. Four turbulence models were explored: SST, Menter’s baseline model (BSL) 
[41], Wilcox’s 1988 k-omega model [55], and a modified version of the SST model by Georgiadis and Yoder 
(SST-GY) [42]. All of these models are variations of the k-omega model with one of the biggest differences 
being the limitation on the turbulent shear stress. The limit is defined as a percentage of the turbulent 
kinetic energy, with 31% for SST, 35.5% for SST-GY, and no limit for BSL and k-omega. For this study, 
BSL and SST-GY models were coded manually by modifying the existing SST model within OVERFLOW. 


0.5 


Ftac Ptete Boundary Conditions, 

M =0,2, Re L = S million (L-lhT* - 540 R 


farf i eld Riemann BC 



. Pt/P-.* 1 . 02028 . 

T t/T ft# = 1.006. 

I quantity from interior 


P/P^hQ — * 
other quantities 
from interior 


symmetry 


adiabanc solid wall 


-0,5 


start of plate at x=0 

i ■ i i * l_ 


1 5 X 1 


L5 


Figure 23: Flat plate 2D grid and boundary conditions [13]. 
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4.2 Results 


A benchmark test with the SST turbulence model was first performed and compared to outputs from CFL3D 
and FUN3D, courtesy of reference [13]. The skin friction coefficient profile for the SST turbulence model for 
all three flow solvers is shown in Fig. 24. It can be seen that the OVERFLOW solution agrees well with the 
other codes. With this confidence, the other turbulence models were run and compared to the SST solution, 
shown in Fig. 24. Experimental data was from Wieghardt and Tillman [56]. It can be seen that the BSL 
and SST-GY models predict a slightly higher skin friction coefficient compared to SST (within 2%) while 
the k-omega model predicts a larger skin friction coefficient compared to SST (within 6%). All models, with 
the exception of k-omega, agree well with the experimental data. 



(a) SST turbulence model. 



Figure 24: Skin friction coefficient profiles. 


To further examine the differences, Fig. 25 shows u + verses y + at a Reynolds number of 4.2 million 
(x = 0.84 ft). It can be seen that all the models agree well with the experimental data and Spalding’s 
formula [57], although the k-omega model under-predicts u+ in the outer boundary layer. This is consistent 
with the over-prediction of the skin fiction coefficient because the upper bound on u+ is inversely proportional 
to the square root of the skin friction coefficient. Note, the following form of Spalding’s formula was used 
[58]: 


y+ = u + + e" 5 - 033 * 


— 1 — kv/ 


1 

2 


(/VM + ) 2 — 


1 

6 


(/VM + ) 


(23) 
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Figure 25: u + verses y + for the various turbulence models at Re x = 4.2E6. 

The stark prediction differences between the k-omega model compared to SST, SST-GY, and BSL are 
more clearly shown in Fig. 26, which includes profiles of the normalized turbulent shear stress and turbulent 
viscosity at the same Reynolds number of 4.2 million. Experimental data in this case was from Klebanoff [59]. 
With the exception of the k-omega model, the CFD solutions agree reasonably well with the experimental 
data. Because of the non-conformity of the k-omega model to the other models tested, an additional study 
varying the freestream turbulent viscosity was performed. The freestream turbulent viscosity was singled 
out because it has been widely known that the 1988 version of the k-omega model is extremely dependent 
on it for a given turbulent kinetic energy state [60]. In fact, one of the motivations behind Mentor’s BSL 
and SST models was to eliminate the freestream turbulent viscosity dependence [41]. The normalized shear 
stress and turbulent viscosity profiles are shown in Fig. 27 from this additional study. Note that the value of 
(Hr/ n) In f used for all prior analyses was 0.3. It can be seen that the k-omega solutions approach the SST 
solution with decreasing freestream turbulent viscosity. However, it most likely would take an unrealisticly 
low freestream turbulent viscosity value for the k-omega model to be within the ballpark of the SST solution 
at a more realistic value. Thus it is not recommended that the 1988 k-omega model be used for the SBLI 
analyses. 
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Figure 26: Normalized turbulent parameters for various turbulence models at Re x = 4.2 E6 




Figure 27: Normalized turbulent parameters at various freestream turbulent viscosities at Re x = A.2E6 
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5 Quasi-ID Code 


5.1 Overview 

For the cases presented, the static temperature of the air in the test section is quite cold at nearly 118 K due 
in part to the high Mach number. Such cold temperatures put the perfect gas assumption in question. In 
particular, the assumption that the specific heat at constant pressure is constant and that the ideal gas law 
holds true. This sensitivity could not be explored by using the production version of OVERFLOW as the 
code cannot handle varying specific heats. Thus, a quasi-ID MATLAB [61] code was developed to compute 
various thermodynamic and flow parameters for calorically perfect and non-ideal air. The latter half required 
the use of the National Institute of Standards and Technology (NIST) Reference Fluid Thermodynamic and 
Transport Properties Database (REFPROP) code [62]. The MATLAB code has the capability to calculate 
the properties at a single station and over an oblique shock for either a given u velocity or A /A* ratio. Total 
conditions must also be provided for the code to run. Station nomenclature in the code uses 1 for upstream 
of the oblique shock and 2 for downstream of the oblique shock. The code is currently set up for air and 
utilizes a predefined mixture [63, 64] containing, by mass fraction, 75.57% nitrogen, 23.16% oxygen, and 
1.2691% argon when interfacing with REFPROP. The code, however, can be modified for a variety of gases, 
such as carbon dioxide. Regardless if the gas is perfect or non-ideal, the code assumes the total enthalpy is 
conserved and that the normal v velocity is constant across the oblique shock, as shown in Figs. 28 and 29. 
The code also assumes single phase states. 



Entropy (kJ/kg-K) 


Figure 28: Sample enthalpy verses entropy diagram for calculating across an oblique shock for a non-ideal 
gas. 1: total conditions at Station 1, IS: static conditions at Station 1, 2S: static conditions at Station 2, 2: 
total conditions at Station 2. 
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5.2 Calorically Perfect Code 

The calorically perfect side of the code (runType options 0 and 10) utilizes the perfect gas equations along 
with the ideal gas law [65]. Station 1 conditions are calculated as follows: 


7 = 1.4 


(24) 


%>, l — 


jR 

7-1 


(25) 


ht,i — 


(26) 


7 7 U 1 

hi = l - y 


(27) 


_ /ii 

Ti = — — 

Cp,i 


(28) 


Mi = 




vWT 


(29) 


Pl = Pt,l ( 1 + 


(30) 
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( 31 ) 


P 1 

Pl RT\ 

To calculate the shock angle and respective Mach number and velocity components ahead of the oblique 
shock [ 66 ]: 


0 = tan 3 (3 + Ci tan /3 + C 2 tan /? + C 3 


(32) 


where 


Cl = C 3 (1 - Ml) 


(33) 


C 2 



Ml + 


2 

7-1 


(34) 


and 


C 3 = 


2 

7- 1 

tan 9 ( Mf + ^x) 


(35) 


u n ^i=uismf3 (36) 

M n? i = Mi sin /3 (37) 

To calculate Station 2 conditions, mass, momentum, and energy equations were balanced across the oblique 
shock. 


Pi^n,! = p2U n ,2 (38) 


Pl + Pl u n,l ~ P2 + P2^n,2 (39) 

2 ^n,l + = 2^,2 + ^2 (40) 
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P 2 = P2RT2 


(41) 


Cp , 2 — %>,1 


(42) 


h 2 = Cp, 2 T 2 


(43) 


Equations (38) through (43) were solved numerically using a bisection scheme. Once converged, the remaining 
Station 2 conditions could be calculated. 


u 2 


u n , 2 


sin (P - 0) 


M 2 


u 2 


V 1 RT 2 


M n 2 = 


2 

ViRT 2 


T t , 2 = T t ,i 


/it, 2 = +,i 


(44) 


(45) 


(46) 


(47) 


(48) 


P/,2 — P2 ( 1 + “7 — M ‘2 


7— 1 A "> 1 

2~ 


(49) 


Once station conditions were found, sonic properties could be calculated for either station. 


A* \7 + 1 


2 \ 2 ( 7-0 


7+1 


1 / 7 — 1 9 \ 2(7_1 ) 

M _1 1 + 2 M 2 


(50) 


r* = t, 1 


7- 1 


(51) 


h* = c p T* 


(52) 
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p* =pt 


7-1 

2 


7 -l 


(53) 


u* = y/jRT* (54) 

7 / 

M* = — (55) 

7X* 

5.3 Non-Ideal Code 

The non-ideal side of the code (runType options 2 and 12) utilizes REFPROP for most of the calculations 
and does not take into consideration the ideal gas law. Station 1 conditions are found as follows: 


h,i~f(T t ,i,pt,i) (56) 

«t,i = f (T t ,i,Pt,i) (57) 

V 2 

hi = h t -^ (58) 

Ti=/(fci,st,i) (59) 

Pi = f (hi,s t ,i) (60) 

/'i = / (7' l l) (61) 


«i = f(hi,s t ,i) 


(62) 


Mi = 


mi 

a x 


(63) 
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71 = / (/' l l ) 


(64) 


c P ,i = f (hi,s t ,i) (65) 

To calculate the shock angle and respective Mach number and velocity components ahead of the oblique 
shock: 


tan (/3 - 0) _ u n ,2 
tan/3 u n , i 

u n ,i = ui sin /3 
M n? i = Mi sin/3 

Equations (66) through (68) were solved numerically using a bisection scheme as an outer-loop to Equations 
(69) through (72), which numerically balance the mass, momentum, and energy equations across the oblique 
shock. 


( 66 ) 

(67) 

( 68 ) 


P l^n,l = P2U n ,2 


(69) 


Pl+PlU 2 n,l =P2 + P2ul t2 


(70) 


1 .2 , u _ 1.2 


-« n> l + /ll = -« n>2 + /l 2 


(71) 


P2 = f{h 2 ,P2) 


(72) 


Once Equations (66) through (72) were converged, the remaining Station 2 conditions could be calculated. 


S2 = f(h 2 ,p 2 ) 


(73) 


T 2 = f{p 2 ,S 2 ) 


(74) 
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a 2 = f(p 2 ,s 2 ) 


(75) 


72 = f(P2,S 2 ) 


Cp, 2 = f (P2,S 2 ) 


— 2 
_ sin (/? - 0) 


Mo = - 


M w o = 


ht, 2 — ^t,i 


Tt.2 = f{h t ,2,s 2 ) 


Pt , 2 — / (^t,2, S 2 ) (83) 

Once station conditions were found, sonic properties could be calculated for either station. Sonic conditions 
were found by first numerically solving for the sonic enthalpy and velocity. 


/i* — rif — 

2 


u* = / (h*, s) 


Once converged, rest of the sonic properties could be calculated. 


P* = f (h*, s ) 
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T* = f (h*, s) 


(87) 


, ^ u 
M* = — 

u* 


( 88 ) 


7* = f(h*,s ) 


(89) 


A p* T 1 
A* p T* M* 


(90) 


5.4 Quasi-ID Exploration 

Several parameters of interest were calculated for a range of Mach numbers, spanning Mach 1.4 to 3.5. The 
percentage difference between these parameters for the calorically perfect and non-ideal air are shown in Fig. 
30. The percentage difference is defined as: 


%Diff = lOOx 


Nonideal — Perfect 
Perfect 


(91) 


Also, the perfect dynamic viscosity values were obtained using the following form of Sutherland’s law [67]: 




1.458xl0 -6 


2 ^ 1.5 

T + 110.4 


(92) 


It can be seen that there is not a drastic difference between the solutions, with at most a 0.4% difference 
between non-ideal and perfect air at Mach 2.75. The only exception to this would be the Prandtl number, 
which is upwards of 9% difference. For this study the perfect air Prandtl number was assumed to be a 
constant value of 0.702, which is associated with air at standard atmospheric conditions. It is known that 
the Prandtl number varies with temperature and therefore most of the difference between the non-ideal and 
perfect Prandtl numbers is due to temperature variance and not from assuming perfect or non-ideal air. 
Also explored was the effect of the shock angle for varying wedge angles, shown in Fig. 31. It can be seen 
that there is not much of a difference between the non-ideal and perfect shock angles with a difference of 
about 0.2% at Mach 2.75, regardless of the wedge angle. The noise within Fig. 31 is most likely due to 
the convergence tolerance when solving Equations (66) through (72). Based on these ID studies it can be 
concluded that although there is some effect from assuming the air is perfect, it is a good enough assumption, 
especially given the extra computational power needed to compute non-ideal air in 3D simulations. 
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% Differance Between Non-Ideal and Perfect Air 



% Differance Between Non- Ideal and Perfect Air 



(b) Zoomed in to d= 0.5%. 


Figure 30: Percentage differences between non-ideal and perfect air for various parameters. 


Shock Angle % Differance Between Non-Ideal and Perfect Air for Varying Wedge Angle 



Figure 31: Percentage differences of the shock angle between non-ideal and perfect air for varying wedge 
angles. 
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6 Grid Resolution Study 

A grid resolution study was performed on the standard SBLI grid with the Standard Case conditions. The 
study looked at solutions obtained on the coarse, medium, and fine grid levels in a similar manor as Galbraith 
[3], with streamwise cross-sections of the grids shown in Fig. 32. Due to OVERFLOW interpolating final 
solutions onto the finest mesh, regardless of grid level, solutions were read into Tecplot [67, 68] by reading 
every fourth point for the coarsest grid level and every other point for the medium grid level. Convergence 
of the solutions was based on the L2 residual, shown in Fig. 33. In addition, the fine grid level solution was 
run out an additional 1000 iterations and showed that the u and v velocity profiles of interest changed by 
no more than 0.1 m/s. Thus, all solutions presented are considered iteratively converged. 


L 



(a) Coarse grid. 

Z 

L. 



(b) Medium grid. 

I 

L. 



(c) Fine grid. 

Figure 32: Streamwise cross-sections at center-span for the various grid levels. 
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(c) Fine grid level. 


Figure 33: L2 residuals for the Standard case. 


Several parameters were looked at to determine grid convergence. First, total pressure, total temperature, 
and Mach number were computed at three different locations, shown in Table 2. Location coordinates can 
be found in Table 4. It can be seen that the totals ahead of the oblique shock are converging to the input 
totals of 98,000 Pa and 295.7 K and that the Mach numbers at the throat and upstream of the oblique shock 
are converging to their ID expected values of 1.00 and 2.75, respectively. Also, it can be seen that the total 
pressure and Mach number downstream of the oblique shock are correctly converging to lower values than 
their upstream counterparts while the total temperature is mostly conserved. 

Second, mass flow rates at the inlet, exit, and throat regions were computed and shown in Table 3. 
Mass flow rates were computed using the mass flow integration routine within Tecplot [67]. Like the total 
quantities and Mach number, the mass flow rates are converged well with less than 0.5% loss throughout 
the tunnel (with 86% of this total mass loss upstream of the throat) on the finest grid level. The mass loss 
is primarily due to round-off error within the discretization process. Finally, u velocity profiles of interest 
within the streamwise plane located underneath the wedge were examined, as shown in Fig. 34. It can be 
seen that while the u velocity profiles are very close for the medium and fine grid levels, they themselves are 
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(wui) a 


not fully grid converged. Because the profiles are crucial for the results presented in this thesis, it is fair to 
say that this parameter is the dominating one out of the parameters presented. Thus while the results are 
shown to be still converging, the fine grid level is considered most suitable for the studies presented here. 




Figure 34: u velocity profiles within the streamwise plane. 


Table 2: Various flow parameters used to determine grid convergence. 



Coarse 

Medium 

Fine 

Pt (Pa) 
Throat T t (K) 

Mach 

98761.8 

296.530 

0.952165 

98199.4 

295.909 

0.947367 

98009.4 

295.701 

0.947031 

Pt (Pa) 

Upstream T t (K) 

Mach 

98134.4 

296.504 

2.73365 

98103.0 

295.960 

2.73637 

97996.9 

295.704 

2.74482 

Pt (Pa) 

Downstream T t (K) 

Mach 

96083.7 

296.173 

2.49111 

96459.3 

296.073 

2.48308 

96247.1 

295.742 

2.47317 


Table 3: Mass flow rates (kg/s) used to determine grid convergence. 



Coarse 

Medium 

Fine 

Inlet 

0.250120 

0.247012 

0.245282 

Throat 

0.246334 

0.245170 

0.244326 

Exit 

0.247058 

0.245823 

0.244169 
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Table 4: Coordinates for grid resolution study. 



X (mm) 

Y (mm) 

Z (mm) 

Throat 

-538.6808 

8.7612 

-2.1299 

Upstream 

-188.1651 

36.4712 

-2.1299 

Downstream 

34.7961 

21.8123 

-0.6699 

Inlet 

-971.9310 

All 

All 

Throat 

-538.6808 

All 

All 

Exit 

403.6530 

All 

All 
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7 Results 


7.1 Standard Case 

The streamwise data plane from the SBLI experiment was used as the primary plane to compare experimental 
and CFD data as it is considered to be the most accurate of all the experimental data planes [19]. The location 
of the plane is shown in Fig. 35, which is located at center-span. Velocity contours for the Standard case 
along with the experimental data are shown in Fig. 36. For consistency with the prior workshop, velocities 
were normalized by Ui n f — 603 m/s. Although the contours appear similar, the CFD solution under predicts 
the velocities (particularly u). In addition, the reflected oblique shock in the CFD solution is upstream of 
the experimental data. To further inspect the flow field, Fig. 37 shows contours of negative u velocities just 
above the bottom wall underneath the wedge. These contours provide a convenient way to approximate the 
bottom wall separation region underneath the wedge. The resulting blockage, derived from the corner flows, 
is typical of prior CFD analysis and has been shown by Galbraith et. al. [23] to play a major role in the flow 
field for this tunnel due to the tunnel’s small size. The blockage at the throat is also important and can be 
quantified by introducing a throat blockage parameter, 6*, which is defined by the following equations [23]: 

m* = (1 - 5*) m* idea i (93) 


where 


rh*ideal = A * 


Pt 


7 


7 + 1 

7—1 


(94) 


VT t \ R V7 + + 

The blockage parameter calculation assumes adiabatic, ideal perfect gas, and choked flow and is shown to 
be 1.53% for the Standard Case. 



X (mm) 

Figure 35: Data plane location in reference to the grid. 
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Experiment 



20 25 30 35 40 45 50 55 60 

X (mm) 


Standard Case CFD 



20 25 30 35 40 45 50 55 60 

X (mm) 


Experiment 



20 25 30 35 40 45 50 55 60 

X (mm) 


Standard Case CFD 



20 25 30 35 40 45 50 55 60 

X (mm) 


Figure 36: u velocity (left) and v velocity (right) streamwise contours for the experiment (top) and Standard 
case (bottom). 



X (mm) 

Figure 37: Bottom wall separation underneath the wedge for the Standard case. 


7.2 Particle Lag Effects 

Velocity contours with the particle lag incorporated are shown in Fig. 38, with the remaining velocity 
contour plots located in appendix A. Note only the Standard and Combined cases are shown as they are 
representative of the spectrum of presented cases. Figures 39 through 42 show u and v velocity profiles with 
and without particle lag at the most upstream axial location along with the first intersecting spanwise data 
plane. See appendix B for remaining axial stations within the streamwise plane. The lag is shown to bring 
both the Standard and Combined cases closer to the experimental data with increasing lag time constant. 
Because the lag has such an effect, it is not advisable to use the workshop metric for comparing CFD to 
this PIV obtained experimental data without accounting for the lag. Recall that it is a fraction of the entire 
experimentally derived time constant that is used for the lag model. Although it is a very short time, it 
creates an error in the velocity field that is higher than the desired match between the CFD and experiment. 
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Y (mm) Y (mm) Y (mm) Y (mm) 


Experiment 


Experiment 


15 

10 

5 



V/Ulnl 

0.1 


- 0.02 

-0.06 

-0.1 


Standard Case CFD, with Short Lag 

15 


V/UInf 



20 25 30 35 40 45 50 55 60 

X (mm) 


Standard Case CFD, with Medium Lag 



20 25 30 35 40 45 50 55 60 

X (mm) 


Standard Case CFD, with Long Lag 
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E 

E 


15 


10 


5 



V/UInf 
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- 0.02 

- 0.06 

-0.1 


Combined Case CFD, with Short Lag 



20 25 30 35 40 45 50 55 60 

X (mm) 


Combined Case CFD, with Medium Lag 



20 25 30 35 40 45 50 55 60 

X (mm) 


Combined Case CFD, with Long Lag 



20 25 30 35 40 45 50 55 60 

X (mm) 


Figure 38: v velocity streamwise contours for the experiment (top), short lag (middle-top), medium lag 
(middle-bottom), and long lag (bottom). 
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Figure 39: Select velocity profiles at the most-upstream axial location (x = 18.191 mm) for the Standard 
case. 




Figure 40: Select velocity profiles at the first intersecting spanwise plane ( x = 20.76 mm) for the Standard 
case. 
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Figure 41: Select velocity profiles at the most-upstream axial location (x = 18.191 mm) for the Combined 
case. 




Figure 42: Select velocity profiles at the first intersecting spanwise plane ( x = 20.76 mm) for the Combined 
case. 


7.3 Isothermal Case 

The effects of the Isothermal case can be seen in Figs. 43 and 44. Figure 43 and the static temperature 
profile in Fig. 44 show little difference in the thermal boundary layer between the Standard and Isothermal 
cases. However, assuming the aluminum is isothermal is more realistic as there is likely heat transfer between 
the aluminum and the tunnel’s surroundings. The temperature difference is evident in the total temperature 
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center-span profile at the bottom wall, shown in Fig. 45. For reference, several notable stations have 
been marked to orient the profile. Figure 45 also confirms that the isothermal boundary conditions were 
implemented correctly and that the deviation from the constant total temperature for the Isothermal and 
Combined cases is due to the presence of the adiabatically modeled bottom window. Further, the isothermal 
boundary conditions were shown to have an effect on the total temperature profile, as shown in Fig. 44. The 
thermal recovery factor, r c , is determined to be 92.3% from the following equation: 


^ = l + §(7-l )Ml (95) 

-Loo ^ 

where T aw is the adiabatic wall temperature, or the total temperature evaluated from an adiabatic calculation. 
This value of r c is reasonable for turbulent flows as the recovery factor is approximately equal to the cube 
root of the Prandtl number. 



(b) Isothermal case CFD. 


Figure 43: Freestream thermal boundary-layer at x = —63.9 mm. Temperature cut off is at 99% freestream 
(121.2 K). 
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20 r 



15 - 



Figure 44: Temperature profiles at the most-upstream axial location {x = 18.191 mm). 



Figure 45: Total temperature centerspan profile at the bottom wall with the following stations marked: 1. 
Throat, 2. Trip location, 3. Start of the tunnel straight section, 4. Wedge leading edge, 5. Wedge trailing 
edge. 


The differences in temperature do correspond to differences in velocity, as shown in Fig. 46. The 
Isothermal case is shown to be slightly worse than the Standard case as compared to the experimental data. 
The differences in the velocities can be seen more clearly in the difference contour plots shown in Figs. 47 
and 48. The difference is defined as the Isothermal case minus the Standard case. It is shown that the 
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Y (mm) Y (mm) 


interaction region for the Isothermal case is further upstream compared to the Standard case, and therefore 
further away compared to the experimental data. The movement of the interaction region is also verified in 
the static temperature difference contour, shown in Fig. 49. 




Figure 46: Select velocity profiles at the most-upstream axial location (x = 18.191 mm). 



20 25 30 35 40 45 50 55 60 

X (mm) 


(a) Positive difference (m/s). 



20 25 30 35 40 45 50 55 60 

X (mm) 


(b) Negative difference (m/s). 


Figure 47: u velocity difference contours. 



X (mm) X (mm) 

(a) Positive difference (m/s). (b) Negative difference (m/s). 

Figure 48: v velocity difference contours. 
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X (mm) X (mm) 

(a) Positive difference (K). (b) Negative difference (K). 

Figure 49: Static temperature difference contours. 


7.4 Trip Case 

The effects of the Trip case are shown in Figs. 50 through 52. Figures 50 and 51 show that Trip and Combined 
cases agree well with the experimental data for the u velocity which suggests that assuming the flow is laminar 
at some portion upstream of the wedge is likely more correct than assuming it to be turbulent throughout. 
This further supports likely transition downstream of the throat based on the momentum thickness Reynolds 
number of the first laminar case. The closeness of the u velocity profiles in Figs. 50 and 51 for the Trip 
and Combined case indicates that the trip is a dominating factor that sets the Combined case apart from 
the Standard case. This is also evident in the turbulent kinetic energy profiles shown in Fig. 52. Having a 
laminar region upstream of the wedge allows for thinner boundary layers, and thus less blockage, upstream 
of the wedge. This is further quantified by the throat blockage parameter, with the Trip case having a throat 
blockage of 1.11%, which is less than the 1.53% that it was for the Standard case. 




Figure 50: Select velocity profiles at the most-upstream axial location (x = 18.191 mm). 
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Figure 51: Select velocity profiles at the first intersecting spanwise plane (x = 20.76 mm). 




Figure 52: Turbulent kinetic energy profiles. 


Examining the v velocity profiles in Fig. 51 shows that the experimental data from the spanwise and 
streamwise planes do not agree with each other. In fact, the CFD solutions are shown to match the experi- 
mental data from the spanwise plane better than compared to the streamwise plane. This is a key point when 
developing an error metric as the error metric used in the prior workshop focused solely on the experimental 
data in the streamwise plane and not the spanwise plane. 
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7.5 Turbulence Modeling Effects 

Figures 53 and 54 show velocity profiles while Fig. 55 shows the turbulent kinetic energy profile for the TKE 
and MUT case. It can be seen that the velocity profiles for the TKE case compare better to the experimental 
data than the Standard case does, however, the TKE case still predicts lower freestream velocities. It can 
also be shown that the dominating factor is the increase in freestream turbulent kinetic energy and not 
the increase in freestream turbulent viscosity. However this domination is in part due to that the inlet 
freestream turbulent viscosity was increased by only lOx while the inlet freestream turbulent kinetic energy 
was increased by 10,000x. The same conclusions can also be made with the turbulent kinetic energy profiles 
of Fig. 55. 




Figure 53: Select velocity profiles at the most-upstream axial location (x = 18.191 mm). 
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Figure 54: Select velocity profiles at the first intersecting spanwise plane (x = 20.76 mm). 




Figure 55: Turbulent kinetic energy profiles. 


Figures 56 and 57 show the velocity profiles for the various turbulence models tested for the Standard 
case conditions. Additional velocity profiles can be found in appendix B. It can be seen that both BSL and 
SST-GY predict higher u velocities in the boundary layer and thus tend to agree better with the experimental 
values compared to the SST case. This is also true for the v velocities within the boundary layer region. 
However, all three turbulence models are the same for the freestream velocities. The greatest impact of the 
turbulence models can be seen in the blockage effects. Figure 58 shows the bottom wall separation underneath 
the wedge. It can clearly be seen that the corner flow separation greatly decreases with turbulence models that 
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allow for greater shear stress transport. However, the center separation bubble is shown to increase slightly 
in the spanwise direction with decreasing corner flow separation. This is consistent with the experimentally 
shown coupling mentioned earlier. It is also shown that the corresponding throat blockage also decreases, 
from the 1.53% with SST to 1.51% with SST-GY to 1.49% with BSL. A similar turbulence model study was 
also performed for the Combined case conditions, which can be found in appendix C. 




Figure 56: Velocity profiles at the most-upstream axial location (x = 18.191 mm) for various turbulence 
models, Standard Case. 




Figure 57: Velocity profiles at the first intersecting spanwise plane ( x = 20.76 mm) for various turbulence 
models, Standard Case. 
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(b) BSL (positive span) and SST (negative span). 


Figure 58: Bottom wall separation underneath the wedge for the Standard case with various turbulence 
models. 


7.6 Metrics 

To compare the CFD data with the prior CFD data from the workshop, the following error metric was used 
[19]: 


e (f)n ~ | {fcfd) n (fexp) n | (96) 

£(/) = V^> (/) « ( 97 ) 

The error metric was summed over the entire streamwise data plane with metric results shown in Table 5. 
Cases submitted to the workshop are designated by letters while the cases presented in this paper are in bold 
face. It can be seen that the Standard case ranks best out of the cases presented and the addition of the lag 
improves both the Standard and Combined cases. Although the metric is a useful quantity for comparison, 
it must be used in combination with all the data. All improvements to matching the data must be for the 
right reasons. Although the Standard case with the PIV lag corrections have a smaller v velocity error, the 
physics of the Combined case is likely more correct. 
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Table 5: Workshop error metric comparisons. 


U Error 

V Error 

0.02373 

Q 

0.008947 

Experiment 

0.02633 

B 

0.01158 

Standard (Medium Lag) 

0.02669 

P 

0.01185 

Standard (Long Lag) 

0.02676 

Standard (Long Lag) 

0.01224 

Standard (Short Lag) 

0.02747 

Standard (Medium Lag) 

0.01308 

MUT 

0.02759 

G 

0.01331 

TKE 

0.02840 

F 

0.01348 

Combined (Medium Lag) 

0.02853 

Standard (Short Lag) 

0.01360 

Combined (Short Lag) 

0.02899 

M 

0.01375 

Combined (Long Lag) 

0.02957 

I 

0.01377 

Standard 

0.02964 

Standard 

0.01403 

Combined 

0.02999 

K 

0.01414 

Trip 

0.03020 

Standard (SST-GY) 

0.01449 

B 

0.03025 

Combined (Short Lag) 

0.01514 

Isothermal 

0.03035 

N 

0.01621 

Modified Geometry 

0.03036 

Combined (Medium Lag) 

0.01682 

P 

0.03043 

TKE 

0.01716 

G 

0.03043 

MUT 

0.01729 

F 

0.03047 

Combined 

0.01771 

M 

0.03064 

Combined (Long Lag) 

0.01828 

Q 

0.03090 

Isothermal 

0.01867 

K 

0.03114 

Modified Geometry 

0.01917 

N 

0.03115 

Standard (BSL) 

0.01950 

Standard (SST-GY) 

0.03129 

0 

0.01961 

O 

0.03163 

Trip 

0.02227 

Standard (BSL) 

0.03473 

Experiment 

0.02344 

j 

0.03571 

H 

0.02348 

Combined (SST-GY) 

0.03739 

Combined (SST-GY) 

0.02576 

Combined (BSL) 

0.03856 

Combined (BSL) 

0.02721 

H 

0.03980 

L 

0.03883 

L 

0.03995 

j 

0.04002 

I 


One possible metric to be used in conjunction with the workshop metric is a comparison of the u velocities 
at points of interest. For this study, two points were chosen: points A and B in Fig. 59. Point A was chosen 
outside the data plane, upstream of the wedge and at center height and span to capture a freestream point. 
Point B was chosen within the data plane at u max at the upstream most axial location. Velocity and 
velocity differences are shown in Table 6. The velocity difference is defined as the current case minus the 
Standard case. Interestingly, the Combined case difference at Point A is reflected as a buildup of the previous 
differences. This is not true at Point B, but as seen in Fig. 59, the flow at Point B is not truly at freestream 
conditions with an apparent decrease in Mach number (and thus decreases in u velocity). This in part 
explains why the workshop error metric predicts the Standard case as being better than the Combined case 
as the data used in workshop error metric is derived from at /downstream of Point B. Also, the highest u 
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velocity at Point B from all the simulations is still 11.6 m/s lower than the experimental value and should 
be just upstream of the shock. 


A: Taken at center- B: Point of most upstream 

height, center-span U max within experimental data 


plane 



Figure 59: Point locations within the centerspan plane. Mach number contour of the Standard case CFD 
solution is shown. 


Table 6: u velocity and velocity differences. 



Point A 

Point B 


U (m/s) 

AU (m/s) 

U (m/s) 

AU (m/s) 

Standard 

594.600 

0.000 

587.042 

0.000 

Isothermal 

594.454 

-0.146 

586.377 

-0.665 

Modified Geometry 

596.567 

1.967 

586.413 

-0.629 

Trip 

595.186 

0.586 

587.737 

0.695 

Combined 

596.980 

2.380 

587.038 

-0.004 

Experiment 

- 

- 

599.330 

12.288 


Another proposed complementary metric would be the use of the oblique shock angle off the leading 
edge of the wedge. The shock angle was calculated by finding the angle between the leading edge of the 
wedge and a point downstream with the max second derivative of a scaled arc length. For details on how 
the shock angle was calculated, see appendix D. The calculated shock angles for the various cases as well as 
for the experimental data are shown in Table 7. It can be seen that the shock angle does not change much 
between the cases, but the addition of the lag brings it to within the experimental value. A comparison to 
quasi-ID theory is shown in Fig. 60. It can be seen that the shock angle cannot be calculated based on the 
quasi-ID theory due to the 3D nature of the shocks present in the flow. The advantage of this metric is that 
separate experimental evaluations of the shock angle are possible: from detailed PIV to optical methods 
such as schlieren. 
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Figure 60: Shock angle verses Mach number. 


Table 7: Shock angle comparisons. 



Shock Angle (deg) 

Standard 

29.888 

Standard (BSL) 

29.888 

Standard (SST-GY) 

29.888 

Isothermal 

29.888 

Modified Geometry 

29.543 

Trip 

29.888 

Combined 

30.231 

Combined (BSL) 

30.231 

Combined (SST-GY) 

30.231 

TKE 

29.888 

MUT 

29.888 

Standard (Short Lag) 

29.196 

Standard (Medium Lag) 

29.196 

Standard (Long Lag) 

28.494 

Combined (Short Lag) 

29.543 

Combined (Medium Lag) 

29.543 

Combined (Long Lag) 

28.846 

Experiment 

29.196 
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7.7 Blockage Revisited 

Figures 61 and 62 show the bottom wall separations underneath the wedge for several of the cases presented. 
Separations are shown to be similar between cases, but there are noticeable differences. In particular, the 
Modified Geometry case was shown to have the most corner flow separation while the Trip case was shown 
to have the least corner separation out of all the cases run with SST. This is further quantified by the throat 
blockage parameter as the Modified Geometry case had 1.56% throat blockage compared to 1.11% of the 
Trip case. The separation regions can also help explain why velocities at Point B (located underneath the 
wedge) tended to decrease compared to the Standard Case while the velocities at Point A (upstream of the 
wedge) increased. The respective increase or decrease is most noticeable for cases with more separation. The 
extent of the separations was thought to be due mostly to details within the turbulence model, however, 
these results show that there are many sources that impact the separations and respective blockages. 



X (mm) 

(a) Isothermal case (positive span) and Standard case (negative span). 



X (mm) 

(b) Modified Geometry case (positive span) and Standard case (negative span). 

Figure 61: Bottom wall separation underneath the wedge for the Isothermal and Modified Geometry cases. 
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(b) Combined case (positive span) and Standard case (negative span). 

Figure 62: Bottom wall separation underneath the wedge for the Trip and Combined cases. 
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8 Conclusions/Future Work 


CFD analyses were performed on the UM Mach 2.75 Glass Tunnel with a 7.75 degree semi-span wedge (a 
test case used at a recent SBLI workshop) to explore key physics pertaining to SBLFs. It was shown that 
the CFD generally under predicted the freestream velocities, however the solutions showed improvement, 
particularly the v velocity component, with several sensitivity insights. These included the addition of a 
laminar region upstream of the wedge with a specified trip point as well as geometry that reflects the current 
state of the tunnel. The flow was shown to be most likely transitional downstream of the throat. This 
is based on the analysis of the momentum thickness Reynolds number for the laminar case as well as the 
improved velocity profiles in the boundary layer upstream of the shock for cases that accounted for the 
transition. Other sensitivities explored included using an isothermal boundary condition, throat geometry 
modifications, and several turbulence models and turbulence model parameter sensitivities. Also, a fraction 
of the measured PIV lag was used with a simple model to modify the CFD solutions and showed improvement 
to the comparison to the PIV results. It showed that the velocity of the PIV particle is not the velocity of 
the air behind the shock, and the difference is larger than one would like for a validation test case. Future 
comparisons should have the CFD results augmented in a post-processing step to calculate particle velocities. 

A variety of metrics were used to capture the results presented in this thesis. One metric used was the 
primary metric from the workshop: a root mean square error norm of the difference between the measured 
u and v velocities. Almost all cases presented in this thesis show improvement in that metric for v velocities 
as compared to the cases run for the workshop. The u velocity errors were better than the experimental 
uncertainty except for two cases. Two new error metrics have been proposed that are felt to be complimentary 
and help establish confidence that the physics is correct. These include a velocity upstream of the interaction 
and the center-span shock angle. 

Although several sensitivities were presented in this thesis, there are still several sensitivities that have 
not been addressed. Additional geometric parameters of interest include the wedge leading edge location and 
angle, wedge width and location within the tunnel, and noticeable wall indents from micro ramp accommo- 
dations. These parameters would shed more insight into the sensitivity of the geometry on the oblique shock. 
Sensitivities to the turbulence model and freestream turbulent kinetic energy should be explored further to 
address the impact of the corner flows, and therefore replicate the correct blockage on the freestream flow 
field. Further explorations of the heat transfer boundary conditions should be explored, such as conjugate 
heat transfer and the sensitivity to changes in wall temperature which is known to change throughout exper- 
imental test runs on a given day. Also, future studies should explore boundary-layer transition/trip location 
as the flow is most certainly transitional downstream of the throat. Aside from the geometric and flow sensi- 


NAS A/TM— 20 13-218081 


56 



tivities, the simplified PIV models should be improved on in order to further understand the effects of PIV as 
related to post-processing CFD data for experimental data comparison. Improvements include calculating 
the particle lag based on the forces exerted on the individual particles by the air (including particle size 
distribution) and obtaining flow field snapshots at two instances in time. These snapshots would then be 
processed using the same PIV post-processing algorithm used with the experimental data. 
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A Velocity Contours 

This appendix contains the u and v velocity contours for the CFD cases presented. 
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Figure A.l: u velocity contours for the Standard (left) and Combined (right) 


cases with lag. 
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Figure A. 2: u (left) and v (right) velocity contours for Isothermal and Modified Geometry cases. 
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Figure A. 3: u (left) and v (right) velocity contours for Trip and Combined cases. 
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Figure A. 4: u (left) and v (right) velocity contours for TKE and MUT cases. 
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Figure A. 5: u (left) and v (right) velocity contours for the Standard case with SST-GY and BSL. 
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Figure A. 6: u (left) and v (right) velocity contours for the Combined case with SST-GY and BSL. 
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B Velocity Profiles 

This appendix contains velocity profiles for all remaining spanwise intersecting planes within the streamwise 
plane. 


B.l Standard Case 




Figure B.l: Velocity profiles at x = 26.76 mm, Standard case. 




Figure B.2: Velocity profiles at x = 30.76 mm, Standard case. 
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Figure B.3: Velocity profiles at x = 34.76 mm, Standard case. 




Figure B.4: Velocity profiles at x = 38.76 mm, Standard case. 
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Figure B.5: Velocity profiles at x = 41.76 mm, Standard case. 




Figure B.6: Velocity profiles at x = 53.76 mm, Standard case. 
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B.2 Combined Case 




Figure B.7: Velocity profiles at x = 26.76 mm, Combined case. 




Figure B.8: Velocity profiles at x = 30.76 mm, Combined case. 
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Figure B.9: Velocity profiles at x = 34.76 mm, Combined case. 




Figure B.10: Velocity profiles at x = 38.76 mm, Combined case. 
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Figure B.ll: Velocity profiles at x = 41.76 mm, Combined case. 




Figure B.12: Velocity profiles at x = 53.76 mm, Combined case. 
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B.3 Modified Geometry, Isothermal, and Trip Cases 




Figure B.13: Velocity profiles at x = 26.76 mm, various cases. 




Figure B.14: Velocity profiles at x = 30.76 mm, various cases. 
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Figure B.15: Velocity profiles at x = 34.76 mm, various cases. 




Figure B.16: Velocity profiles at x = 38.76 mm, various cases. 
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Figure B.17: Velocity profiles at x = 41.76 mm, various cases. 




Figure B.18: Velocity profiles at x = 53.76 mm, various cases. 
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B.4 TKE and MUT Cases 




Figure B.19: Velocity profiles at x = 26.76 mm, TKE and MUT cases. 




Figure B.20: Velocity profiles at x = 30.76 mm, TKE and MUT cases. 
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Figure B.21: Velocity profiles at x = 34.76 mm, TKE and MUT cases. 




Figure B.22: Velocity profiles at x = 38.76 mm, TKE and MUT cases. 
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Figure B.23: Velocity profiles at x = 41.76 mm, TKE and MUT cases. 




Figure B.24: Velocity profiles at x = 53.76 mm, TKE and MUT cases. 
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B.5 Standard Case with Various Turbulence Models 


See appendix C for Combined case equivilant. 




Figure B.25: Velocity profiles at x — 26.76 mm, Standard case with various turbulence models. 




Figure B.26: Velocity profiles at x = 30.76 mm, Standard case with various turbulence models. 
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Figure B.27: Velocity profiles at x = 34.76 mm, Standard case with various turbulence models. 




Figure B.28: Velocity profiles at x = 38.76 mm, Standard case with various turbulence models. 
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Figure B.29: Velocity profiles at x = 41.76 mm, Standard case with various turbulence models. 




Figure B.30: Velocity profiles at x = 53.76 mm, Standard case with various turbulence models. 
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C Combined Case with Various Turbulence Models 


This appendix contains velocity profiles and bottom wall separation plots for the Combined case with the 
SST-GY and BSL turbulence models. Similar trends were noticed compared to the turbulence model study 
with the Standard case conditions. However unlike the Standard conditions, the throat blockage was inconsis- 
tent with the turbulent shear stress limiter, with blockage values of 0.930% with SST, 0.997% with SST-GY, 
and 0.988% with BSL. This is most likely because unlike the Standard case conditions, the Combined case 
conditions modeled the throat region laminar and thus should be independent of turbulence quantities. The 
discrepancy between the turbulence models shown is therefore due to numerical error. 




Figure C.l: Velocity profiles at x = 18.191 mm, Combined case with various turbulence models. 
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Figure C.2: Velocity profiles at x = 20.76 mm, Combined case with various turbulence models. 




Figure C.3: Velocity profiles at x = 26.76 mm, Combined case with various turbulence models. 
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Figure C.4: Velocity profiles at x = 30.76 mm, Combined case with various turbulence models. 




Figure C.5: Velocity profiles at x = 34.76 mm, Combined case with various turbulence models. 
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Figure C.6: Velocity profiles at x = 38.76 mm, Combined case with various turbulence models. 




Figure C.7: Velocity profiles at x = 41.76 mm, Combined case with various turbulence models. 
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Figure C.8: Velocity profiles at x = 53.76 mm, Combined case with various turbulence models. 



X (mm) 

(a) SST-GY (positive span) and SST (negative span). 



X (mm) 

(b) BSL (positive span) and SST (negative span). 


Figure C.9: Bottom wall separation underneath the wedge for the Combined case with various turbulence 
models. 
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D Shock Angle Calculation 

This appendix outlines how the shock angle was computed for the various cases. One would think computing 
a shock angle would be a trivial exercise, but that is not the case. The method outlined below utilizes v 
velocity profiles since they vary dramatically after being turned through a wedge angle. The method adopted 
is then consistent for the experimental as well as the CFD results. The shock angle was defined as the angle 
of the line formed by two points, denoted as the shock line, and shown in Fig. D.l. The start point of the 
shock line was chosen to be the wedge leading edge in the streamwise plane as this point is fixed for all 
cases examined. To find the end point of the shock line, a v velocity profile was first obtained at an axial 
station of interest. In this case, the axial location was chosen to be the most upstream location within the 
streamwise data plane (x = 18.191 mm). For consistency with the experimental data, the v velocity profile 
was obtained by interpolating the CFD data onto the PIV grid. Second, the v velocity profile was scaled by 
the profile aspect ratio to obtain the profile quantity, x: 


, Umax Umin 
X = V 


Vinax Vmin 


Third, the arc length profile was computed, with the arc length defined as: 


(D.l) 


h = 


\J (z'i +1 - x 'if + (y i+ i - yif 


+ h-1 


(D.2) 


Forth, the first and second derivatives of the arc length with respect to y were computed using a forward 
differentiating scheme. 


f dl\ li -\- 1 li 

\dy)i Vi+i-yi 


(D.3) 



( 4l\ _ f dl_\ 

\ d y)i + 1 \ d y)i 

Vi+i ~ Vi 


(D.4) 


7 2 7 

The point with the maximum ^ became the end point for the shock line. The shock angle was then 
computed as: 


(tf) (D - 6) 

For a visual example, see Fig. D.2 through Fig. D.4, which show the above mentioned profiles for the 
Standard case. 
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Figure D.l: Shock angle diagram/nomenclature. 
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(b) x profile at x = 18.191 mm. 


Figure D.2: Velocity profiles for determining shock line end location (Standard case shown). 


35 

30 


+ CFD 
+ Max(d2l/dy2) 


25 


20 

15 



0 24 6 8 10 12 

Y (mm) 


16 


18 


Figure D.3: Arc length profile for determining shock line end location (Standard case shown). 
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Figure D.4: Derivative profiles for determining shock line end location (Standard case shown). 
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E Simulation Checklists 


Checklists for the simulations were used to document each run in a consistent manner as outlined by Galbraith 
[3]. 

E.l Initial 


Input File 

Spatial order (FSO) 

1 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

4 

CFL 

0.5-2 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

SST 

Turbulent Input 
Quantities 


K 

IE-2 

Epsilon 

0.1 

Reference 

Quantities 


Mach 

2.75 

Tref (R) 

205.771 

Re 

636160 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

None 

Check reference 

/ 

conditions 

RefMach 

2.75 

XKInf 

IE-2 

Retlnf 

0.3 

RefL (m) 

0.06985 

Ttot (K) 

292.778 

Ptot (Pa) 

101125 

Check the initial 

/ 

solution 

Remove orphan 
nodes 

/ 

Use most resent 

/ 

INTOUT file 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

2000, 2000, 2000 

Dimensionalize the 
solution 

/ 

Comments 

Unconverged 
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E.2 Standard Case 


Input File 

Spatial order (FSO) 

3 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

1 

CFL 

5-20 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

SST 

Turbulent Input 
Quantities 


K 

IE-6 

Epsilon 

0.3 

Reference 

Quantities 


Mach 

2.75 

Tref (R) 

211.845 

Re 

592877 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

Initial 

Check reference 

/ 

conditions 

RefMach 

2.75 

XKInf 

IE-6 

Retlnf 

0.3 

RefL (m) 

0.06985 

Ttot (K) 

295.7 

Ptot (Pa) 

98000 

Check the initial 

/ 

solution 

Remove orphan 
nodes 

/ 

Use most resent 

/ 

INTOUT file 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

0, 0, 20000 

Dimensionalize the 
solution 

/ 

Comments 

Converged 
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E.3 Isothermal Case 


Input File 

Spatial order (FSO) 

3 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

1 

CFL 

5-20 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

SST 

Turbulent Input 
Quantities 


K 

IE-6 

Epsilon 

0.3 

Reference 

Quantities 


Mach 

2.75 

Tref (R) 

211.845 

Re 

592877 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

Standard Case 

Check reference 

/ 

conditions 

RefMach 

2.75 

XKInf 

IE-6 

Retlnf 

0.3 

RefL (m) 

0.06985 

Ttot (K) 

295.7 

Ptot (Pa) 

98000 

Check the initial 

/ 

solution 

Remove orphan 
nodes 

/ 

Use most resent 

/ 

INTOUT file 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

0, 0, 20000 

Dimensionalize the 
solution 

/ 

Comments 

Converged 
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E.4 Modified Geometry Case 


Input File 

Spatial order (FSO) 

3 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

1 

CFL 

5-20 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

SST 

Turbulent Input 
Quantities 


K 

IE-6 

Epsilon 

0.3 

Reference 

Quantities 


Mach 

2.75 

Tref (R) 

211.845 

Re 

592877 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

Initial 

Check reference 

/ 

conditions 

RefMach 

2.75 

XKInf 

IE-6 

Retlnf 

0.3 

RefL (m) 

0.06985 

Ttot (K) 

295.7 

Ptot (Pa) 

98000 

Check the initial 

/ 

solution 

Remove orphan 
nodes 

/ 

Use most resent 

/ 

INTOUT file 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

0, 0, 20000 

Dimensionalize the 
solution 

/ 

Comments 

Converged 
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E.5 All-Laminar Case 


Input File 

Spatial order (FSO) 

3 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

1 

CFL 

5-20 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

SST 

Turbulent Input 
Quantities 


K 

IE-6 

Epsilon 

0.3 

Reference 

Quantities 


Mach 

2.75 

Tref (R) 

211.845 

Re 

592877 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

Initial 

Check reference 

/ 

conditions 

RefMach 

2.75 

XKInf 

IE-6 

Retlnf 

0.3 

RefL (m) 

0.06985 

Ttot (K) 

295.7 

Ptot (Pa) 

98000 

Check the initial 

/ 

solution 

Remove orphan 
nodes 

/ 

Use most resent 

/ 

INTOUT file 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

0, 0, 20000 

Dimensionalize the 
solution 

/ 

Comments 

Converged 
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E.6 Trip Case 


Input File 

Spatial order (FSO) 

3 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

1 

CFL 

5-20 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

SST 

Turbulent Input 
Quantities 


K 

IE-6 

Epsilon 

0.3 

Reference 

Quantities 


Mach 

2.75 

Tref (R) 

211.845 

Re 

592877 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

All-Laminar Case 

Check reference 

/ 

conditions 

RefMach 

2.75 

XKInf 

IE-6 

Retlnf 

0.3 

RefL (m) 

0.06985 

Ttot (K) 

295.7 

Ptot (Pa) 

98000 

Check the initial 

/ 

solution 

Remove orphan 
nodes 

/ 

Use most resent 

/ 

INTOUT file 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

0, 0, 10000 

Dimensionalize the 
solution 

/ 

Comments 

Converged 
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E.7 Combined Case 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

Modified Geometry 
Case 

Check reference 
conditions 

/ 

RefMach 

2.75 

XKInf 

IE-6 

Retlnf 

0.3 

RefL (m) 

0.06985 

Ttot (K) 

295.7 

Ptot (Pa) 

98000 

Check the initial 
solution 

/ 

Remove orphan 
nodes 

/ 

Use most resent 
INTOUT file 

/ 


Input File 

Spatial order (FSO) 

3 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

1 

CFL 

5-20 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

SST 

Turbulent Input 
Quantities 


K 

IE-6 

Epsilon 

0.3 

Reference 

Quantities 


Mach 

2.75 

Tref (R) 

211.845 

Re 

592877 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

0, 0, 20000 

Dimensionalize the 
solution 

/ 

Comments 

Converged 
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E.8 TKE Case 


Input File 

Spatial order (FSO) 

3 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

1 

CFL 

5-20 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

SST 

Turbulent Input 
Quantities 


K 

IE-2 

Epsilon 

0.3 

Reference 

Quantities 


Mach 

2.75 

Tref (R) 

211.845 

Re 

592877 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

Initial 

Check reference 

/ 

conditions 

RefMach 

2.75 

XKInf 

IE-2 

Retlnf 

0.3 

RefL (m) 

0.06985 

Ttot (K) 

295.7 

Ptot (Pa) 

98000 

Check the initial 

/ 

solution 

Remove orphan 
nodes 

/ 

Use most resent 

/ 

INTOUT file 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

0, 0, 20000 

Dimensionalize the 
solution 

/ 

Comments 

Converged 
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E.9 MUT Case 


Input File 

Spatial order (FSO) 

3 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

1 

CFL 

5-20 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

SST 

Turbulent Input 
Quantities 


K 

IE-2 

Epsilon 

3.0 

Reference 

Quantities 


Mach 

2.75 

Tref (R) 

211.845 

Re 

592877 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

Initial 

Check reference 

/ 

conditions 

RefMach 

2.75 

XKInf 

IE-2 

Retlnf 

3.0 

RefL (m) 

0.06985 

Ttot (K) 

295.7 

Ptot (Pa) 

98000 

Check the initial 

/ 

solution 

Remove orphan 
nodes 

/ 

Use most resent 

/ 

INTOUT file 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

0, 0, 20000 

Dimensionalize the 
solution 

/ 

Comments 

Converged 
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E.10 Standard Case with SST-GY 


Input File 

Spatial order (FSO) 

3 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

1 

CFL 

5-20 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

SST-GY 

Turbulent Input 
Quantities 


K 

IE-6 

Epsilon 

0.3 

Reference 

Quantities 


Mach 

2.75 

Tref (R) 

211.845 

Re 

592877 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

Standard Case 

Check reference 

/ 

conditions 

RefMach 

2.75 

XKInf 

IE-6 

Retlnf 

0.3 

RefL (m) 

0.06985 

Ttot (K) 

295.7 

Ptot (Pa) 

98000 

Check the initial 

/ 

solution 

Remove orphan 
nodes 

/ 

Use most resent 

/ 

INTOUT file 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

0, 0, 20000 

Dimensionalize the 
solution 

/ 

Comments 

Converged 
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E.ll Standard Case with BSL 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

Standard Case with 
SST-GY 

Check reference 
conditions 

/ 

RefMach 

2.75 

XKInf 

IE-6 

Retlnf 

0.3 

RefL (m) 

0.06985 

Ttot (K) 

295.7 

Ptot (Pa) 

98000 

Check the initial 
solution 

/ 

Remove orphan 
nodes 

/ 

Use most resent 
INTOUT file 

/ 


Input File 

Spatial order (FSO) 

3 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

1 

CFL 

5-20 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

BSL 

Turbulent Input 
Quantities 


K 

IE-6 

Epsilon 

0.3 

Reference 

Quantities 


Mach 

2.75 

Tref (R) 

211.845 

Re 

592877 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

0, 0, 5000 

Dimensionalize the 
solution 

/ 

Comments 

Converged 
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E.12 Flat Plate with SST 


Input File 

Spatial order (FSO) 

3 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

1 

CFL 

5-40 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

SST 

Turbulent Input 
Quantities 


K 

IE-6 

Epsilon 

0.3 

Reference 

Quantities 


Mach 

0.2 

Tref (R) 

540 

Re 

5E6 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

None 

Check reference 

/ 

conditions 

RefMach 

0.2 

XKInf 

IE-6 

Retlnf 

0.3 

RefL (ft) 

1 

Ttot (K) 

302.39 

Ptot (Pa) 

386192 

Check the initial 

/ 

solution 

Remove orphan 
nodes 

/ 

Use most resent 

/ 

INTOUT file 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

5000, 5000, 220000 

Dimensionalize the 
solution 

/ 

Comments 

Converged 
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E.13 Flat Plate with K-Omega 


Input File 

Spatial order (FSO) 

3 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

1 

CFL 

5-40 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

K-Omega 

Turbulent Input 
Quantities 


K 

IE-6 

Epsilon 

0.3 

Reference 

Quantities 


Mach 

0.2 

Tref (R) 

540 

Re 

5E6 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

Flat Plate with SST 

Check reference 

/ 

conditions 

RefMach 

0.2 

XKInf 

IE-6 

Retlnf 

0.3 

RefL (ft) 

1 

Ttot (K) 

302.39 

Ptot (Pa) 

386192 

Check the initial 

/ 

solution 

Remove orphan 
nodes 

/ 

Use most resent 

/ 

INTOUT file 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

0, 0, 180000 

Dimensionalize the 
solution 

/ 

Comments 

Converged 
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E.14 Flat Plate with SST-GY 


Input File 

Spatial order (FSO) 

3 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

1 

CFL 

5-40 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

SST-GY 

Turbulent Input 
Quantities 


K 

IE-6 

Epsilon 

0.3 

Reference 

Quantities 


Mach 

0.2 

Tref (R) 

540 

Re 

5E6 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

Flat Plate with SST 

Check reference 

/ 

conditions 

RefMach 

0.2 

XKInf 

IE-6 

Retlnf 

0.3 

RefL (ft) 

1 

Ttot (K) 

302.39 

Ptot (Pa) 

386192 

Check the initial 

/ 

solution 

Remove orphan 
nodes 

/ 

Use most resent 

/ 

INTOUT file 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

0, 0, 180000 

Dimensionalize the 
solution 

/ 

Comments 

Converged 
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E.15 Flat Plate with BSL 


Input File 

Spatial order (FSO) 

3 

Time order (TFOSO) 

1 

Time-step scaling 
flag (ITIME) 

1 

CFL 

5-40 

Turbulence CFL 
(CFLT) 

1 

Turbulence order 
(FSOT) 

1 

Turbulence model 

BSL 

Turbulent Input 
Quantities 


K 

IE-6 

Epsilon 

0.3 

Reference 

Quantities 


Mach 

0.2 

Tref (R) 

540 

Re 

5E6 


Pre-Processing 

Orient meshes 

/ 

No negative volumes 

/ 

Check boundary 
conditions for gaps 

/ 

Initial Solution 

Flat Plate with SST 

Check reference 

/ 

conditions 

RefMach 

0.2 

XKInf 

IE-6 

Retlnf 

0.3 

RefL (ft) 

1 

Ttot (K) 

302.39 

Ptot (Pa) 

386192 

Check the initial 

/ 

solution 

Remove orphan 
nodes 

/ 

Use most resent 

/ 

INTOUT file 


Post Processing 

Iterations (Coarse, 
Medium, Fine) 

0, 0, 180000 

Dimensionalize the 
solution 

/ 

Comments 

Converged 
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